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We study weakly first order cosmological phase transitions in finite temperature field theories. 
Focusing on the standard electroweak theory and its minimal supersymmetric extension, we identify 
the regimes of Higgs masses for which the phase transition in these models proceeds by significant 
phase mixing and the coarsening of the subsequent domain network. This dynamics is distinct from 
that for strongly first order transitions, which proceed by the nucleation and propagation of critical 
bubbles. We describe how electroweak baryogenesis might take place in these models, explaining 
how our new picture can relax the sphaleron washout bound of traditional scenarios. 
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Finite temperature phase transitions are of great in- 
terest in cosmology because they provide a mechanism 
by which remnants of the early universe can make an 
observable imprint on today's universe. Well-known ex- 
amples of this are cosmological inflation, the production 
of topological defects, and electroweak baryogenesis . 

Cosmological phase transitions are well-understood in 
two particular limits. Strongly first order phase transi- 
tions occur when a system begins in a metastable state 
that is separated from a global minimum by a sufficiently 
large energy barrier. In such a system, widely separated 
points in space undergo quantum tunneling or thermal 
hopping events, in which bubbles of the true vacuum 
nucleate in the background sea of false vacuum: bub- 
bles which are sufficiently large that their volume energy 
dominates over their surface tension expand and even- 
tually coalesce, completing the phase transition. If the 
phase transition is second order (continuous), the sys- 
tem begins in an unstable state, with no energy barrier 
separation from the global minimum. Small-amplitude, 
large-wavelength fluctuations grow, followed by domain 
coarsening and phase separation. This is referred to as 
spinodal decomposition 

Away from these two extreme limits, the dynamics of 
phase transitions is much less well understood. Numeri- 
cal simulations provide accurate estimates of the strength 
of the phase transition, and determine the point at which 
first order transitions become second order as parameters 
of the theory are varied, as recently done for the standard 
model (SM) However, for weakly first order phase 
transitions an understanding of how the phase transition 
proceeds and completes is far more vague. 

It is widely believed that, in the absence of exotic 
means for departing from equilibrium, a strongly first 
order phase transition is required for electroweak baryo- 
genesis to occur. (For a recent review see |^.) This is 
for two main reasons. If bubble nucleation is the mech- 
anism, then the sharp change in order parameter across 
the bubble walls provides a large departure from equilib- 
rium at each point in space swept out by the walls. Fur- 



ther, the associated large energy barrier between phases 
leads to thermal anomalous fermion number violation be- 
ing sufficiently suppressed in the true vacuum so that 
any baryon number produced is not washed out 1^. For 
a weakly first order transition this is not the case, and 
the traditional criterion that the transition be sufficiently 
strong provides a constraint on the theory. This con- 
straint can be expressed as {4>{Tc)) /Tc > 1, where (f) is 
the Higgs field, and Tc is the critical temperature for 
the electroweak phase transition. This translates directly 
into a bound on the mass of the Higgs boson, which, for 
the minimal supersymmetric standard model (MSSM) is 
^(MSSM)^^Q5 GeV@. 

In this letter, we investigate an alternative scenario 
for fermion production in weakly first order phase tran- 
sitions. When transitions are weak, large-amplitude sub- 
critical fluctuations between the symmetric and broken- 
symmetric phases cause the transition to begin by sig- 
nificant phase mixing and, after inter-phase fiuctuations 
cease, to proceed by coarsening of the subsequent domain 
network. The crucial point is that the sphaleron washout 
temperature need not be the critical temperature, but the 
lower temperature when the large-amplitude fluctuations 
freeze out, the Ginzburg temperature, Tq. Thus, fermion 
production and preservation may still be efficient in this 
picture, and possibly remain sufficient to account for the 
baryon asymmetry of the universe. Let us describe how 
this might happen using the SM and the more promising 
MSSM as examples. 

The finite temperature effective potential for the mag- 
nitude of the Higgs field in electroweak models can often 
be written as 



V{cj), T) = D{T^ - T^)(j)^ - ET(f>^ + ^^(f>^ 



with T| = {mjj - 8Bv'^)/AD, where, for the SM, 
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and A(T) is the temperature-corrected Higgs coupling. It 
is meaningful to speak of there being a global minimum, 
separated from a metastable minimum by an energy bar- 
rier, for T| < T2 < Tl = j^^T^. Let us begin by ex- 
amining how the sphaleron washout temperature varies 
in this range, for changing Higgs mass. In figure (^) we 
plot {(f){T))/T as a function of the Higgs mass mn, for 
three values of the temperature, Tc, Tq, and T2. 

FIG. 1. Standard Model Washout. 
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The lower curve on this figure demonstrates how the 
the traditional sphaleron bound arises, since this curve 
crosses unity at a Higgs mass around 40 GeV. However, 
this, along with the upper curve shows that, for Higgs 
masses in the range 40 GeV < mf/^ < 48 GeV, the 
temperature below which sphaleron washout is avoided 
may be reached before the barrier between phases has 
vanished. This is the relevant regime for our scenario. 

Although the SM is useful for illustrating this effect, 
experimental bounds place the physical Higgs mass at 
77iif > 100 GeV. Thus, our analysis is academic in this 
model. However, in the MSSM, we shall see that useful 
results are obtained. For calculational convenience, in 
this letter we study the MSSM in the limit in which the 
extra Higgs fields and all the sparticles are decoupled 
from the spectrum. Thus, all that remains is the lightest 
Higgs and the right-handed stop. In this limit, the one- 
loop finite-temperature effective potential can be written 
in the same form as the SM (|^), but with the following 
parameters 
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For this potential, we can again examine how the 
sphaleron washout temperature varies with Higgs mass. 
We note that two-loop corrections to the MSSM potential 
have been extensively discussed in the literature. In fact, 
the consensus is that the inclusion of two-loop effects en- 
hances the strength of the transition, further opening the 
window of allowed Higgs masses for efficient baryogenesis 
[0. As our present purpose is to illustrate a new possi- 
ble mechanism for baryogenesis within weak first-order 
phase transitions, we will for simplicity restrict ourselves 
to one-loop corrections here, since there will always be 
a range of Higgs masses for which the transition will be 
weakly first order. However, we are presently working to 
extending our analysis to include two-loop corrections. 
We plot our results for the one-loop MSSM potential in 
figure (I). 



FIG. 2. 
Washout. 
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Notice that, for the MSSM, there is a range of masses 
for the lightest Higgs, within which the phase transition is 
weakly first order, but with the possibility for sphaleron 
transitions to become suppressed in the broken phase be- 
fore the end of the transition. 

The phase transition begins at the critical tempera- 
ture, Tc- At this temperature, when the phase transition 
is sufficiently weakly first order (the sphaleron bound is 
violated at Tc but not at Tq), there are significant sub- 
critical fluctuations, both of the broken-symmetric (the 
true vacuum, and the symmetric one (false vacuum, 
0). Consequently, significant phase mixing occurs and ho- 
mogeneous nucleation breaks down Q , |^ . We write the 
rates of the relevant subcritical fluctuations as Go^+ and 

G-i >o, respectively. As the phase-mixed plasma cools 

due to the expansion of the universe, it reaches the tem- 
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perature Tq, at which fluctuations back to the false vac- 
uum freeze out. The criterion for this to occur is 
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where H is the Hubble parameter. The question of how 
the phase transition proceeds after this depends on how 
the fraction, 7-|.(rG), of the total volume in the bro- 
ken phase at the Ginzburg temperature compares with 
the percolation threshold, Ppcrc — 0.31, in three dimen- 
sions. If 7+(Tg) < Pperc, we refer to partial phase mix- 
ing initially. In this case, isolated domains of the -I- 
phase will grow due to pressure difference. However, if 
7_i_(T(5) > Ppcrc, then we have total phase mixing initially, 
and both phases percolate, being separated by a convo- 
luted wall, plus small isolated domains of each phase. 

For both scenarios, baryon- number production will de- 
pend on the ratio {(I)(Tq)) /Tq , which is shown in fig- 
ures (1) and (2) as the (central) dashed line. Clearly, 
the details of how baryon number will percolate into the 
broken-symmetric phase is determined by how 7+(Tg) 
compares to Ppcrc, that is, on whether the broken- 
symmetric phase has percolated or not. However, our 
focus here is on what we could call the necessary condi- 
tion for baryogenesis to occur, which is fixed by the ratio 
{(^{Tq)) /Tg, irrespective of the details of the phase tran- 
sition. There are two possibilities, described as follows: 
If this ratio is larger than one, baryon-number violation 
is suppressed in the broken phase at freeze-out {Tq) and 
a net baryonic excess will be generated as domains of the 
broken phase advance on the symmetric phase. If this ra- 
tio is smaller than one, any baryon excess will be erased 
in the broken phase until the temperature drops so that 
the ratio {(j){T)) /T reaches unity, say at a temperature 
T. Thus, in this second case we may expect a severe 
(if not deadly) volume-suppression factor, unless the do- 
mains advance sufficiently slowly so that Tq > T: if Tq 
were not very close to T, the domains of the symmetric 
phase - the sources of baryon-number excess, would have 
plenty of time to shrink into oblivion before T could be 
reached. 

Clearly, we require an estimate of Tq and ^+(Tq). To 
proceed, we assume that the fluctuations are spherical, 
and consider fluctuations of amplitude 0oc (</'+c) and spa- 
tial size, R> the correlation length, with profiles 



^+^,(r) = (^+,-^+)e-'-'/«'+, 



(5) 



This ansatz has been shown to be quantitatively in agree- 
ment with numerical simulations ||l0[| . These fluctuations 
are drawn from a Gibbs distribution, 



where -F+ is the free energy of the respective fluctuation. 

The prefactor A is expected to be small, within the range 
1 < A< 100, although its calculation for subcritical bub- 
ble nucleation is still an open question. A given fluctu- 
ation (j) away from a local minimum of the free-energy 
density has a free energy "cost" 
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With the ansatz above, these free energies take the form 
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with Ci = [D(T^ - T|) - 3ET(j)+ + (3/2)A</.2], and C2 
i-ET + X(b+). 

The correlation length in the broken phase, 
[V"i<j>+)]-\is 



^ = AD{T^ - T 
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where = 9E^T^ /[8XD{T^ - r|)] > 1. Thus, with (|), 
(p^), and (pT|), we may calculate Tq via (^. 

What remains is to compute 7+(Tg). We proceed by 
considering the Boltzmann equation for the production 
of bubbles of the two competing phases, a procedure valid 
below percolation. This is given by flOfl: 



^+ ' 3-/f = -\v\:^ + {l- 7+)Go^+ - 7+G+^o 



dt 



dR 
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where v is the shrinking velocity of the subcritical bub- 
bles, a{t) is the cosmic expansion factor, and the volume- 
fraction in the broken phase is given by 
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where = dn+/dRd4> is the distribution function of 
disconnected domains (hence the d) of the broken phase 
modeled by the subcritical bubbles, and n+ is their num- 
ber density. This equation can be integrated by defin- 
ing the quantity Y+ = ff/s, where s = ^^g^^T^, is 
the entropy density, and looking for equilibrium solutions 
Y+ = 0. The solution can be written as 



7+ 



-^1 (0max; -^min) 



max; -^min) ^" ^2(0max: -^min ) 



(14) 



where 
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Using (□)-(^), we get an expression for the volume frac- 
tion 7+(Tg). An analytical solution can be obtained if 
we neglect the cubic term in (^) , which is a good approx- 
imation for small enough bubbles, far from percolation. 
The remaining integrals are then Gaussian, and can be 
performed to give 
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^ and e,„ax = 8r/(3^V2^0^,,) is the 



where, p 

thermal length-scale for the system which emerges from 
the calculation. An identical expression holds for /2, but 
with 0max replaced by A+ ee |0,nax - 0+|, and e,nax 
replaced by 9+^^ = %T / {Z^:^/2^^?A_ The ratio A/\v\ 
can be obtained numerically, as in [ |lO| or, in principle, 
analytically, although we leave it as a free parameter here. 

FIG. 3. 7+(rG) for the standard model. 
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Using these expressions, we can estimate 7+(Tg) for 
the SM and for the MSSM. We plot the results in figures 
(^) and (||), respectively. 

We can now refer back to figures (l)-(4) and analyze 
the different scenarios. There are two different issues that 
should be discussed: first, the mass range for which the 
condition (0(Tg))/Tg > 1 is satisfied; and second, how 
7(Tg) compares with the percolation threshold proba- 
bility, Pperc = 0.31. For the MSSM, the case of most 
immediate interest, {4>{Tq))/Tq > 1 for m 
GeV (Figure 2), while j{Tg) 
GeV, for A/\v\ = 400 (Figure 4). Thus, in this scenario, 
the broken-symmetric phase does not percolate within 
the range of Higgs masses for which {(piTc)) /Tq > 1. 
The allowed mass range for baryogenesis is extended to 
at least m^^^^ = 116 GeV. This is a lower bound be- 
cause, as we remarked earlier, even if {4>{Tg)) /Tq < 1, it 
is possible for the condition {(j){T))/T > 1 to be achieved 
at a lower temperature T, as long as T is fairly close to 
Tg (otherwise the transition would be completed before 
T and no net baryon excess would be produced). Thus, 
it is possible, in principle, to extend the allowed mass 
range even further than the value dictated by the con- 
dition {(I){Tg)) /Tg > 1, although a detailed analysis of 
the dynamics of the phase transition is required in order 
to make a more quantitative statement. An analogous 
calculation for the SM yields an allowed mass range ex- 
tended to at least = 45 GeV. We should comment 
on the choices for the values of ^/|ti|. The values for the 
prefactor of the subcritical bubble nucleation rate {A) 
or its shrinking velocity in a plasma (v) are not known. 
However, from results for critical bubble nucleation and 
expansion in the electroweak transition, it is reasonable 
to expect that 1 < ^ < 100 and 0.005 < t; < 0.5. We 
are treating ^/|f| as a free parameter (fortunately, only 
the ratio of the two is relevant), and values between 100 
and 1000 seem to be reasonable, judging from the ranges 
of A and Hopefully, detailed simulations of weakly 
first order transitions could further restrict the range of 
A/\v\. 

We have described how electroweak baryogenesis can 
proceed even if the phase transition is weakly first order. 
This idea is radically different to existing lore about such 
scenarios. In weak phase transitions, baryon number vi- 
olating processes are sufficiently copious in the broken 
phase that any baryon excess is washed out at Tc- How- 
ever, the important point is that the transition dynamics 
can be such that, at a lower temperature, when baryon 
number violation becomes exponentially suppressed in 
the broken phase (at T), the transition may not yet have 
completed. There then remains the possibility for a suffi- 
cient baryon excess to be produced during the remainder 
of the transition. Thus, this approach may have impor- 
tant ramifications for electroweak baryogenesis. In this 
letter we have described how to estimate the relevant 
quantities in these scenarios analytically. In a later pub- 
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lication we will present a detailed semi-analytical and 
numerical approach in which we will obtain more precise 
values. This will involve taking into account the full form 
of the free energy, the effects of nonabclian gauge inter- 
actions, and an accurate analysis of the phase transition 
dynamics and freeze-out rates. 

FIG. 4. 7+(Tg) for the MSSM. The two curves are for 
A/\v\ = 400 (continuous line) and 1000 (dotted Une). 
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